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Abstract 

The leading face of Saturn's moon lapetus, Cassini Regio, has an albedo only one 
tenth that on its trailing side. The origin of this enigmatic dichotomy has been 
debated for over forty years, but with new data, a clearer picture is emerging. 
Motivated by Cassini radar and imaging observations, we investigate Soter's 
model of dark exogenous dust striking an originally brighter lapetus by modeling 
the dynamics of the dark dust from the ring of the exterior retrograde satellite 
Phoebe under the relevant perturbations. In particular, we study the particles' 
probabilities of striking lapetus, as well as their expected spatial distribution 
on the lapetian surface. We find that, of the long-lived particles (> 5 /xm). 
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most particle sizes 10 /im) are virtually certain to strike lapetus, and their 
calculated distribution on the surface matches up well with Cassini Regio's 
extent in its longitudinal span. The satellite's polar regions are observed to 
be bright, presumably because ice is deposited there. Thus, in the latitudinal 
direction we estimate polar dust deposition rates to help constrain models o f 



2010 ). 



thermal migration invoked to explain the bright poles (jSpencer and Denkl h 
We also analyze dust originating from other irregular outer moons, determining 
that a significant fraction of that material will eventually coat lapetus — perhaps 
explaining why the spectrum of lapetus' dark material differs somewhat from 
that of Phoebe. Finally we track the dust particles that do not strike lapetus, 
and find that most land on Titan, with a smaller fraction hitting Hyperion. As 
has been previously conjectured, such exogenous dust, coupled with Hyperion's 
chaotic rotation, could produce Hyperion's roughly isotropic, moderate-albedo 
surface. 

Keywords: lapetus,, irregular satellites, Saturn, satellites, debris disks, 
celestial mechanics 



2 



1. Introduction 



Over three dozen dark irregular satellites have been discovered around Sat 



urn using ground-based telescopes ( 


Gladman et al. 


2001 


Sheooard et al. 


2003 


Jewitt et al. 


2005 




Sheooard et al. 


2006 


). Numerical simulations show that 



these irregular satellites must have undergone intense coUisional evolution that 



would have generated 



system (|Nesvornv et al 



large 



2003 



quantities of dark d ust over the age of th e solar 



Turrini et al 



20091 ). Indeed, 



Bottke et al. 



poioh 



estimate that on the order of 10^*^ kg of dust a thousandth the mass of the 
Earth's moon) has been generated in the outer Saturnian system t hrough coUi- 



siona l grinding of these satellites. Furthermore, the recent discovery ([Verbiscer et al. 



20091) of a vast dust ring originating from the largest of the irregulars (Phoebe) 



shows that these dust-producing collisional processes are ongoing even today. 



Small dust particles are strongly affected by radiation forces ([Burns et al 



19791) : in particular, Poynting- Robertson drag will cause particles to lose energy 
and slowly migrate toward their parent planet. One should therefore expect 
mass transfer from the dark outer irregular satellites to the generally brighter 
inner regular satellites (see Fig. [T]). lapetus is the out ermost of the regula r 



19711) 



satellites and, importantly, is observed to be tidally locked (iMcCord et al 
As such, one hemisphere permanently faces the direction of motion and would 
plow through the cloud of dark dust as the cloud evolves inward. Since Phoebe 
(and most of the other irregulars) orbits retrograde and would generate dust 
particles on retrograde paths, collisions with the prograde lapetus would occur 
at high relative velocities (~ 7 km/s) and the dust would mostly coat lapetus 
only on its leading side. This mod el for t h e exog enous origin of the dark material 
on lapetus was first proposed by ISoteii (|l974l ) and seems plausible; when one 



looks at the observed albedo map of lapetus, one finds that the dark region, 
Cassini Regio, is centered precisely around the apex of motion — a difficult fact 
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Figure 1 : This schematic diagram depicts the expected extent of the Phoebe ring (solid gray) , 
as well as the orbits of both Phoebe (black) and the rest of the irregular satellites (light 
gray). The circle at the center represents the main rings, with the E ring surrounding them . 
Approximate scale is provided in Saturn radii (Rs). According to the model of Sote^ ||1974| ). 
as the ring of retrograde dust drifts inward from Phoebe's semimajor axis (arrows), lapetus' 
leading side sweeps up this material, darkening its leading side. Figure provided by Matthew 
S. Tiscareno. 



to exp lain for endogenou s mechanisms (jPenk et al 



As 



Denk et al 



2010 ) 



()2010l) point out, however, the extremely sharp boundaries 
between bright and dark material cannot be the result of simple dust deposition, 
which would yield more g radual transitions. To exp lain the striking boundaries, 



as well as the bright poles, 



Spencer and DenkI (|201Cll ) propose a model of runaway 



ice sublimation in which areas initially darkened by dust become completely 
blackened. The sublimed ice then settles on the poles and on the brighter 
(and therefore colder) trailing side. These two processes together, exogenous 
dust deposition coupled with thermal ice migration, seem the most promising 
mechanism to forming lapetus' striking global albedo dichotomy. 
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20091) represents a 



The recently discovered "Phoebe Ring" (jVerbiscer et al. 
snapshot in time of the inexorable process of mass transfer from the dark outer 
irregulars onto the inner icy satellites. The ring thickness implicates Phoebe 
as the source, showing that Soter's mechanism of coating lapetus is ongoing. 
We will show below that almost all particles in the Phoebe ring of size > lO/im 
will strike lapetus; however, particles smaller than ^ bfim will strike Saturn, its 
main rings, or esc ape the system within a half-Saturn orbit(~ 15 yrs) due to 



radiation pressure (|Verbiscer et al 



20091) . 



Our presentation improves on two previously published works. 



Burns et al. 



and 



(|l996 ) published a sh ort analysis of the dynamics of dust particles from Phoebe 



Tosi et al. 



(|2010[ ) included in their paper a simplified analysis that included 
Poynting-Robertson drag but neglected the important effects of the dominant 
component of solar radiation pressure (radial from the Sun) , which affects par- 
ticles' eccentricities and can quickly drive small grains out of the system. 

This paper performs a more in-depth analysis, considering all the important 
radiation and tidal perturbations from the Sun and calculating the expected 
coverage on the lapetus surface. We also include the precession of lapetus' 
orbital axis, which extends coverage over the poles. 

The paper is organized as follows. Sec. 2 discusses the determination of prob- 
abilities for dust striking lapetus from numerically integrated dust orbits. In 
Sec. 3 we present calculated distributions of dust on the lapetus surface, com- 
paring them to the observed distribution and using them to obtain estimates 
for polar deposition rates. Sec. 4 addresses the same process for the dozens of 
irregular satellites other than Phoebe, and Sec. 5 tracks the fate of dust grains 
that do not strike lapetus and that instead collide with Hyperion and Titan. 
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2. Collision Probabilities 

2.1. Orbital Integrations for Dust Particles 

In order to estimate dust particles' likelihoods of striking lapetus, we first 
consider the important effects of the perturbations affecting dust particle dy- 
namics. Most dust particles spend their lifetimes in a radial range (between the 
orbits of Phoebe and lapetus) where the dominant perturbations are solar. The 
important modifications to these particle orbits therefore result from radiation 
pressure and solar gravity. Nevertheless, since small particles' eccentricities can 
bring them closer to Saturn, we also included the perturbation from Saturn's 
second-order zonal harmonic in our numerical integrations. 

As mentioned in the introduction, particles smaller than ^ 5/im are so af- 
fected by solar radiation pressure that they are quickly removed from the system; 
in this size regime, electromagnetic forces from the planet's magn etosphere are 



20011) . Note that 



negligible relative to the other perturbations (I Burns et al. 
the smallest particles might not be blown out by radiation pressure; once the 
particle size becomes small relative to the incident light's wavelength, the dust 
particles will no longer be able to effectively couple to the radiation field. Such 
particles presumably account for a small fraction of the total mass, and their 
orbits would decay too slowly to reach lapetus — we therefore ignore them. 
The equation of motion can be written as 



GMs.SAQpr^ SA 
r = ^ — rH —S 



(1) 

where the terms, in sequence, are due to the dominant Saturnian gravity, so- 
lar radiation pressure, Poynting-Robertson drag, the Sun's tidal gravity, and 
Saturn's 32. G is the gravitational constant, Ms Saturn's mass, r the dust 
particle's distance from Saturn, S the solar flux at the particle's position, A the 
particle's cross-sectional ctrccL, Qpr the grRin's pressure efficiency, tti the pa-rticle 
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mass, c the speed of light, a the semi-major axis of Saturn (assumed to be on 
a circular orbit about the Sun), Rs the radius of Saturn, J2 Saturn's second- 
order zonal harmonic, and P2 the second Legendre polynomial. The vector r 
is the particle's velocity, and the other vectors can be seen in Fig. [SJ f is the 
direction from Saturn to the particle's position, S is the direction from the Sun 
to the particle position, a is the direction from the Sun to Saturn, and § is the 
direction along Saturn's spin axis (perpendicular to the equatorial plane). 

Saturn Saturn Spin Axis 

Sun a 



Particle 




Figure 2: Schematic diagram showing the geometry of the important perturbations acting on 
dust grains in orbit around Saturn. Vectors are described in the text above. 



While commonly considered in a heliocentric context, Poynting- Robertson 
drag also causes particle s' orbits around a host planet to decay into the planet 



on a timescale given by (|Burns et al 



19791) : 



TPR ^530 years X (2) 
Pr/g 

where asat is Saturn's semimajor axis in AU (w 9.5) and /S^/g is the dimen- 
sionless ratio of the radiation force to the Sun's gravitational force (in this case 
approximately 0.36/r, where r is the particle size in /im). rp/? therefore scales 
linearly with particle size. 

Superimposed on this slow orbital decay (timescale ^ 10^ years for spherical 
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10 /im particles) is a fast oscillation in the eccentricity (P '--^ 1 Saturn 



due to both solar radiation press ure and the Sun's tidal gravity (IBurns et al 



1979; 



Hamilton and Krivov 



yr c± 30 yrs 



19961 ). Eventually, dust-particle orbits will cross 



that of lapetus as Poynting- Robertson drag reduces the orbit size and radiation 
pressure periodically induces large eccentricities. Over time, therefore, the dark 
particles will impact lapetus' leading side. 

As opposed to gravitational accelerations, accelerations due to radiation 
forces are mass — and therefore size — dependent. As a result, we numerically in- 



tor "dl" (see 



Hamilton 


1993 


Hamilton and Krivov 


1996 


Hamilton and Kriiger 



20081 ). This provides a particle's orbital elements as a function of time for each 
particle size. 

In any particular history, we choose particles of a given size and assign 
them a density (we assume that dust particles would share Phoebe's density 
of 1.6 g/cm?). As discussed in further detail below, they are then started at 
various positions along Phoebe's orbit and initially move with Phoebe's velocity. 
We determine that for our assumed density, particles smaller than 4 /im are 
so affected by radiation pressure that within the first half-Saturn year their 
eccentricities reach a value of unity and the grains either collide with Saturn 
or its rings, or escape the Saturn system entirely. This corresponds to 10 
particle orbits and a negligible probability of collision with lapetus. One should 
therefore expect only a significant contribution to lapetus from particles > 4 
/im in size. Since dust particles are not actually spherical and will contain some 
void space, our assumed density is probably high and our 4 /xm likely represents 
a lower limit. 

On the other extreme, the orbital eccentricities of particles larger than 500 



/xm (Poynting-Robertson decay timescale > 50 million years) are almost com- 
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pletely unaffected by radiation forces and are dominantly affected by the Sun's 
tidal gravitational force, which is independent of particle size. We therefore 
run integrations for particle sizes of 5, 10, 25, 50, 100, and 500 /im. A 25 //m 
particle's orbital element evolution is shown in Fig. [3] with its slow semimajor 
axis decay and rapid eccentricity oscillations(P ^ 30 yrs). The bottom panel 
shows the pericenter distance q. When q crosses a satellite's semimajor axis, 
collisions with that moon become possible. 
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Figure 3: Evolution of a 25/xm particle under the effects of solar perturbations. Top panel 
shows the particle's semimajor axis, which starts at Phoebe (a 2l5Rs) and decays on a 
timescale ~ 2.5 Myrs. Superimposed on this slow evolution of the semimajor axis is a rapid 
oscillation in the eccentricity (middle panel) on a timescale Ri 1 Saturn year fx 30 yrs. The 
bottom panel shows the particle orbit's pericenter q, along with the semimajor axis of lapetus 
and Titan. 



A few considerations supply the appropriate initial conditions for the integra- 
tions. All particles leaving Phoebe must have initial speeds > Phoebe's escape 
speed Vesc- Since dust-producing impact events produce a distribution of cjecta 
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velocities with a decaying tail toward higher speeds, one shou ld expect most 



parti cles that escape Phoebe to have launch speeds near Vesc (jFarinella et al 



19931) . Therefore, since Phoebe's escape velocity is much smaller than its orbital 
velocity (~ 0.1 km/s vs. ^ 1.7 km/s), we expect most dust particles generated 
in an impact with Phoebe to approximately share that moon's orbital elements. 
This sets the initial conditions for the semimajor axis, eccentricity and inclina- 
tion (a = 1.296x10^ km, e — 0.156, i — 175.2° with respect to Saturn's orbital 
plane about the Sun). 

The last three initial conditions — the three angles that determine the orienta- 
tion of the orbit (see Fig. H]) — depend on the time of impact itself. Specifically, 
they are set by the orientation of Phoebe's orbit (51 and w), and Phoebe's po- 
sition within its orbit (/, the true anomaly) at the time of impact. This would 
represent a formidable phase space to cover for long integrations, but fortunately 
several considerations limit the phase space considerably. 

The shortest timescale for the perturbations involved is the ~ 30- yr period 
of the Sun's apparent motion about Saturn. Since the dust particles' orbital 
periods around Saturn are much shorter than this (w 1.5 years), the exact 
position of Phoebe (/) in its orbit at the time of impact does little to infiuence 
the subsequent evolution of the orbit's shape or orientation; thus it can be 
chosen arbitrarily. 

Phoebe's orbit orientation at the time of impact, however, is important 
because it precesses more slowly. Nevertheless, one can still limit the phase space 
by inspecting the geometry. Since the dominant perturbations in this problem 
are all of solar origin, the logical plane from which to reference inclinations is 
Saturn's orbital plane (i.e., the plane in which the Sun appears to move in a 
Saturnocentric frame). Phoebe's inclination relative to Saturn's orbital plane 
of 175° means its own orbital plane is almost coplanar, albeit in a retrograde 
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ORBIT ^ 




Figure 4: The t hree Eulcr angles that define an orbit's orientation: i, Q and w. (Figure from 
lGreenberelll982l') . 



sense, with this reference plane. 

In the limit of coplanarity, only one angle (rather than both D, and cj) is 
required to specify the orientation of the orbit given the inclination, i.e., the 
angle between an arbitrary reference direction and the orbit's pericenter. In 
this case, the physically meaningful reference direction is the one toward the 
source of perturbations, the Sun. The orbital evolution of the dust particle 
therefore does not depend strongly on fl and uj independently, but rather on 
the combination - w, which specifies the angle from the Sun's direction to 
pericenter. Note that for a prograde orbit, the angle from the Sun's direction 
to pericenter would he uj ~ + lu, but since uj is measured in the direction of 
orbital motion, the appropriate combination for retrograde orbits is 51 - 
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The approximations discussed above transform an intractable multidimen- 
sional space of initial conditions into a simple one-dimensional space. The ele- 
ments a, e, and i are those of Phoebe's orbit, and the only other initial condition 
left to supply is the quantity 57 — with fl measured relative to the Sun's di- 
rection at the time of impact. Since impacts could happen at any point in 
the precession cycle, we chose to perform integrations for eight equally-spaced 
values of — w. 

We therefore generate, for each particle size, eight sets of a(t), e(t), and 
i(t) corresponding to eight equally-spaced initial values of 51 — oj. Taking the 
initial values of 51 — w as equally likely, we average over the eight sets of outputs, 
yielding, for each particle size, a single set of functions a(t), eft), and i(t). These 
provide the inputs for the collision probability calculations. While we exploit 
Phoebe's orbit's near-alignment with Saturn's orbital plane to combine 51 and 
uj for our initial conditions, the numerical integrations are carried out fully in 
three dimensions. This allows us to track the orbital inclination, a crucial input 
to a 3-D collision probability calculation. 

2.2. Collision Probabilities 

In or der to estimate c ollision probabilit ies, we used the forma l ism d evel- 



oped by 



Greenbera (|l982l ). as improved by iBottke and GreenbersJ (|l993l ). In 



this formalism, the dust particle's and lapetus' semimajor axes, eccentricities 
and inclinations {a, e, i) are taken as known while the precession angles that 
determine the orientation of the orbits (51 and uj) are treated as uniformly dis- 
tributed. Barring resonances between the orbital periods of the dust particles 
with lapetus, this should be a good assumption over collision timescales (> 10^ 
years), which are long compared to the longest precession timescale (lapetus' 
orbit pole, r ^ 10'^ years). This assumption was found to agree with the angular 
distributions from the numerical integrations. 
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An alternative strategy could have been to numerically integrate many dust 
particle orbits and to directly see when and where on lapetus they strike. One 
drawback of such a method is that lapetus' small size relative to the dust orbits 
would dictate using extremely small step sizes in the integration. Our approach 
allowed us to perform fewer computationally expensive orbit integrations per 
particle size in exchange for computationally cheaper collision probability inte- 
grals. 

The calculations of lGreenbera (jl982l) are too complicated to reproduce here. 



The calculation is performed by first calculating the values of (ilp, fJ/, Wp, cj/) 
that would lead to the two orbits crossing ('p' subscripts refer to the particle 
and T subscripts to lapetus). Then the objects' finite size is taken into ac- 
count by Taylor-expanding around these crossing solutions to find the volume 
in (ilp, il/, cjp, LJi) space over which collisions are possible. Finally one calculates 
from Keplerian theory the probability that both objects will simultaneously be 
close enough to the point of closest approach for a collision to occur within 
one object's orbit. The ratio of this probability to the orbital period provides 
a collision frequency. We c ompared our code to the test cases presented in 



Bottke and Greenbera (|l993r ) and found it reproduced their results. 

Due to the wide disparity between orbital period {'^ 1 year) and the colli- 
sion timescale 10^ years), it is impractical to calculate collision probabilities 
for every orbit. One can see in Fig. [3l however, that while the eccentricity 
is oscillating rapidly, the envelope that bounds the oscillation changes slowly, 
on roughly the Foynting- Robertson timescale {rpji > 1 million years). In par- 
ticular, the figure uses only 10* equally-spaced points in time and is still able 
to capture the full behavior. As a result, rather than calculating probabilities 
every orbit, we did so for 10'* timesteps. 

Given a, e, i for both lapetus and dust particle, Greenberg's formalism 
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(1982) provides a collision frequency, 



Frequency 



Probability of collision within one orbit 



(3) 



Period of orbit 



For timesteps At "C 1 / Frequency^ one can then straightforwardly express the 
collision probability within At as 



One can then recursively generate a cumulative probability of collision C, i.e., 
the probability at time t that the particle has already struck lapetus. Starting 
with C(0) = 0, 



The probability of two collisions within a single At is negligible and was ignored. 
P{ti) depends on the orbital elements for lapetus and the dust particle at ti. For 
the dust particle we used a(ti), e(ti), and i(ti), generated as described above. 
For lapetus, we used the present values of a = 3.561x10® km and e = 0.03. 
lapetus' inclination with respect to Saturn's orbital plane, however, changes 
significantly over time and must be considered more carefully. 

To first approximation, the orbit normal precesses uniformly and at a con- 
stant inclination to a vector determined by the perturbations causing the pre- 
cession. This causes the orbit normal to sweep out a cone (see FiglSJ. The 
vector around which orbits precess defines the local Laplace plane (normal to 
this vector). At Phoebe's orbit, all the dominant perturbations are solar, so the 
local Laplace plane corresponds to the plane in which the Sun appears to move, 
Saturn's orbital plane. Close to Saturn, where the dominant perturbation is 
Saturn's oblateness, the local Laplace plane is Saturn's equatorial plane. lape- 
tus has the unique orbital property among satellites of existing at a distance 
where Saturnian and solar perturbations are comparable, and the local Laplace 



P = Frequency x At. 



(4) 



C{t,) = C(t,_i) + (1 - C(i,_i)) * P{t,). 



(5) 
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plane is intermediate, at about 11.5° to Saturn's orbital plane (see IWardlll98ll ) 



Saturn Orbit Norma! 

A ILPN 



, ION 



Saturn 
Orbital 
Plan 




PON 



Figure 5: A schematic representation of the changing orientations of lapetus' and Phoebe's 
orbits (represented by their respective orbit normals PON = Phoebe Orbit Normal and ION 
= lapetus Orbit Normal). The moons' orbit normals precess at constant inclinations (5° and 
8° for Phoebe and lapetus, respectively) to the normal vector to their local Laplace planes, 
sweeping out a cone. Phoebe's Laplace plane coincides with Saturn's orbit normal, while 
lapetus' local Laplace plane normal (ILPN) is inclined about 11° to Saturn's orbit normal. 



Because of this misalignment, although lapetus will precess at approximately 
constant inclination to the normal to its local Laplace plane, its inclination 
relative to our reference plane (Saturn's orbital plane) will change as the orbit 
precesses (see Fig. O. We therefore assumed uniform precession and coarsely 
averaged the probability calculation over an entire precessional cycle, sampling 
more finely when the inclinations of the particle and lapetus were antiparallel 
and the collision probability was changing fastest. 

Finally, as mentioned in Sec. 2.1, we averaged over the eight equally probable 
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initial conditions that we integrated, yielding an overall cumulative probability 
of collision for the given particle size. 

Apart from lapetus, we also tracked collisions with Hyperion and Titan, 
as well as re-impacts into Phoebe. We therefore straightforwardly generalized 
the discussion above to not only update the cumulative probability of collision 
with lapetus at each timestep, but also those with the other three moons. As 
is discussed below, Titan's large size renders it a sink for any long-lived dust 
particles that cross its path; thus, no other moons interior to it would receive 
appreciable amounts of dust and such bodies are therefore not tracked. As 
mentioned earlier, however, a significant fraction of the particles smaller than 
~ 5 fim whose eccentricities all reach unity will strike Saturn or its rings within 
the first half Saturn-year. 

2.3. Results 

Figure [S] shows the calculated cumulative collision probability with lapetus 
for 5, 10, and 25 /im grains. 

Cumulative Protaatality of Collision with lapetus vs. time 
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Figure 6: Cumulative collision probabilities vs. time for 5, 10 and 25 /im particles. Particles 
> 10/im almost all strike lapetus, though larger particles take a longer time to do so. 

Particles 10 /im and larger, being less affected by radiation forces, evolve 
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inward via Poynting-Robertson drag so slowly (i.e., they execute many lapetus- 
crossing orbits before crossing the orbits of Hyperion or Titan) that they almost 
all eventually strike lapetus. As stated before, particles smaller than about 4 /im 
quickly strike Saturn or escape the Saturn system. Therefore, of the longer-lived 
particles, almost all particle sizes are bound for lapetus — only a very narrow size 
range (between about 4 and 10 /im for the chosen density) can miss and end up 
mostly on Titan, with a substantially smaller fraction striking Hyperion. 

We mention that, while almost all particles larger than ^ 10 /im would 
eventually strike lapetus, it takes larger particles longer to evolve inward and hit 
the satellite. In particular, for particles in the geometrical optics limit (the peak 
wavelength in the solar spectrum ^ 0.5/im <C 2'Krdust, satisfied for all particle 
sizes we consider) , the Poynting-Ro bertson decay timescale grows linearly with 



particle size (see 



Burns et al 



19791) . As a reference, assuming particles share 
Phoebe's density of 1.6 g/cm^, 10 /im particles reach lapetus in w 1 Myr. 

As particle sizes increase, one should expect to find a threshold where parti- 
cles stop hitting lapetus when the Poynting-Robertson decay timescale becomes 
longer than the timescale for the destruction of dust grains. Unfortunately, de- 
str uction lifetimes for dust in the outer Saturn system are not well constrained 
(cf. iBurns et al 



20011 ). One mechanism for the destruction of dust grains is 
through mutual collisions. One can estimate the mean free time between parti- 
cle collisions as 

P 

tMF ^ — (6) 
T 

where P is the particles' orbital period and r the ring's norma l optical depth. 



Taking the optical depth in the Phoebe ring, r 2 x 10 ^ (jVerbiscer et al 



20091 ) ■ this yields Imf 100 Myr, the Poynting-Robertson decay timescale 
corresponding to 1 mm grains; however, for each particle size, only collisions 
with particles of roughly the same size or larger affect the dynamics. This 
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would act to increase tMF, but is dependent on the (currently unconstrained) 
particle size distribution. On the other hand, dust rings co Uisionally generate d 



20101) 



early in the Solar System likely had higher optical depths (jBottke et al. 
lowering tMF- For this work we chose the maximum upper-size cutoff imposed 
by setting the Poynting-Robertson decay timescale equal to the lifetime of the 
Solar System. This yields a particle size of ~ 1 cm. Improved estimates of 
coUisional dust lifetimes in the Phoebe ring must await further observations. 

As the introduction mentions, a large supply of dust has been available in 
the outer Saturn system over the course of the Solar System's history. The fact 
that particles > 10 /im are virtually certain to strike lapetus strongly implicates 
coUisionally generated dust as the trigger to lapetus' stark albedo dichotomy. 
An exogenous origin of the dark material explains why the pattern is centered 
on the apex of lapetus' motion and, as shall be shown in Sec. 3, the dynamics 
predict a wrapping of dark material onto the trailing side consistent with that 
observed. 

2.4- Titan, the gatekeeper to the inner Saturnian system 

We now explain the sharp drop in the final fraction of particles that strike 
lapetus between 5 and 10 /im, as seen in Fig. [51 This is due to the large ec- 
centricities induced by radiation pressure, visible in Fig. [S] For the smallest 
particles, the eccentricities are high enough that before the dust grains' prob- 
abilities of striking lapetus near certainty, their orbits begin to cross that of 
Titan. Saturn's largest moon is such a better interceptor of particles that the 
probability of striking lapetus quickly stops increasing and levels off. 

There are several reasons why Titan is highly efficient at eliminating dust 
particles. Most obviously, its sheer size makes its geometrical cross section 
larger than lapetus' by a factor of about 12. Another reason is that collision 
rates depend on the objects' relative velocity, as this determines how frequently 
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the objects can potentially encounter each other (see discussion following Eq. 
I32p . Relative velocities between dust particles and Titan are substantially higher 
than those with lapetus simply because in order to reach the further-in Titan, 
particles generally have to be on very eccentric orbits (e — 0.7-0.9), and will 
encounter Titan close to periapse. 

One might have expected slow relative velocities to lead to enhanced colli- 
sion probabilities due to strong gravitational focusing for slow encounters. In 
fact, gravitational focusing plays little role in this problem because the moon 
orbits are prograde while the dust orbits are retrograde, resulting in high rela- 
tive velocities compared to the satellites' escape velocities (wesc — 0.572 km/s 
for lapetus, v^sc — 2.639 km/s for Titan). For typical encounter velocities, 
lapetus' gravitational cross-section is about 0.5% greater than its geometrical 
cross-section (< 10% for Titan). 

Though we account for gravitational enhancements to the collision cross 
section, in our orbit integrations we ignore close encounters with Titan (and all 
other satellites) on subsequent orbital paths. For typical relative velocities, the 
maximum scattering angle from a close encounter is w 10°. The corresponding 
angle for lapetus is ~ 0.5°. 

We postpone our discussion of the total amount of material that strikes 
Titan and the smaller Hyperion, along with its implications, until Sec. 5. 

3. Coverage 

Since much of the dust previously orbiting in the outer Saturnian system 
will eventually strike lapetus, we ask where on lapetus those particles would 
have landed. In particular, can the dynamics match the extent of lapetus' dark 
side, Cassini Regio? The emplacement of dust could then trigger the thermal 
migration of ice thought to give lapetus the striking appearance it has today 



(jSpencer and Denk 



20101) . 
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Cassini Regio extends beyond lapetus' leading sid e by tens of d egree s onto 



Burns et al 



(|l996l ) have 



the trailing side along the equator (see Fig. [7]). As 
suggested, dust eccentricities naturally explain the longitudinal extension of the 
dark material onto the trailing side. 

If orbits were perfectly circular, dust particles would only strike lapetus 
head-on, as both objects would be moving perfectly azimuthally in opposite 
directions; thus, only the leading face would be darkened since, as discussed 
at the end of Sec. 2.4, the encounter velocities make gravitational focusing 
negligible. 

When particles have eccentric orbits, however, particle velocities are no 
longer perfectly azimuthal, and the radial components allow particles to strike 
the moon further along the equator (see Fig. [S]). Eccentricities induced by 
radiation pressure therefore provide a natural mechanism for extending dust 
coverage onto the trailing side. 
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Figure 7: Global mosaic of lapetus (from lAlberil2008l '). Dark Cassini Regio is centered around 
the apex of lapetus' motion, roughly at 90°W, and extends tens of degrees beyond 0° and 
180° W onto the trailing side. The bright poles (beyond ~ ±60° latitude) an d sharp bo undaries 
between light and dark terrain are likely the result of thermal ice migration l lSpencer and Dcnd 
l2010t) . 



However, just as eccentricities act to extend coverage longitudinally, dust- 
orbit inclinations and lapetus' varying orbital tilt should extend coverage lati- 
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Figure 8: lapctus is depicted as the circle moving on a prograde orbit, while the dust moves 
on retrograde orbits. The white lines separate lapetus' leading and trailing sides. When 
orbits are circular (left), dust will solely darken the leading side, while the radial velocities of 
eccentric orbits allow dark material to reach part of the trailing side. 



tudinally over the poles (j Burns et al 



19961 and see Fig. [5|). Images of lapetus, 



however, reveal bright, icy poles. 

As previously men tioned, thermal ice migr ation provides a mechanism for 



brightening the poles ([Spencer and Denk 



2010|). Icy patches on lapetus dark- 



ened by exogenous dust increase in temperature as a result of their lowered 
albedo. Sublimation rates, which depend exponentially on temperature (e.g.. 



Vvazovkin and Wight 



19971) , thereby increase sharply. This liberates bright ice 
and leaves behind an even darker surface. The further darkened surface's tem- 
perature rises further, and the cycle repeats in a self-accelerati ng process until a 



lag de posit forms with thickness of order the thermal skin depth 



Spencer and Denk 



(|201Cl ). The result is that warm, darkened areas become extremely dark and 
ice-free, while the sublimed ice settles on the coldest areas of the moon — the 
trailing side and the poles. 



The distribution of dark material on the surface therefore holds several in- 
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sights into ongoing processes on lapetus as well as to the past and present 
prevalence of dark dust in the outer Saturnian system. Unfortunately, it is diffi- 
cult to observationally determine the dark layer's depth. Bright-floored craters 
from small impactors that punctured through the dark layer constrain th e layer 



to being much thinne r than the crater's depth ^ 10 m (IDenk et al 



2010!) , while 



20061 ) imply Cassini Regio is on the order of 



radar measurements (lOstro et al. 
decimeters deep. 

Given the above background, we wish to calculate the probability distribu- 
tion for where on lapetus dust would strike for three reasons: 

a) One can convert a probability distribution to a depth distribution (Sec. 
3.2) and compare the resulting global map to the observed lapetus surface. Such 
a comparison tests the hypothesis that lapetus is darkened by dust from Phoebe 
and can provide depths in areas where observations are not available. 

b) Calculated polar deposition rates of dust yield an estimate of the mini- 
mum sublimation rate required to overwhelm dust deposition and keep the poles 
bright. 

c) A global depth distribution provides the total volume of dark material on 
lapetus. This volume, coupled with the collision probabilities of dust calculated 
in Sec. 2, provides a probe of the total amount o f dust collisionally generated 



Bottke et al 



2010|). 



in the outer Saturnian system over its history (cf. 

We subdivide this problem by first calculating the collision probability distri- 
bution over the surface of lapetus in Sec. 3.1. Then 3.2 converts this probability 
distribution to a depth distribution, and 3.3 estimates the sublimation rates re- 
quired to keep the poles bright. We postpone discussion of point c) to Sec. 
5. 
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3.1. Collision Probabilities as a Function of Latitude and Longitude 

We now find the probability per unit area for particles striking lapetus at 
latitude 9 and longitude (/). Note that we can quickly determine the rough 
shape such a distribution should take in the limit of circular, uninclined orbits 
(a good approximation for large particles). In this limit, dust particles strike 
lapetus' leading side head-on. Also, since a/^i?/, we can approximate the 
orbits of lapetus-striking particles as parallel straight lines. Finally, in this 
approximation, our assumed uniform distribution in the variables D, and uj for 
both orbits (see Sec. 2.2) translates into a uniformly distributed bundle of 
quasi-parallel trajectories capable of striking lapetus. In such a uniform field, 
the probability of an impact in a given area element simply is proportional to its 
projected area, given by dAcosip, where ip is the angle between lapetus' velocity 
vector and the outward normal vector to the area element. Equivalently, is 
the angular distance from the apex of motion (see Fig. [9l in which lapetus is 
moving to the left). In this simple case then, the probability per unit area is a 
simple function of ip, 

P{e, 4>) cx cos ip. (7) 

This approximation is good over most of the leading hemisphere, though it 
is clearly incapable of describing the extension of the dark material onto the 
trailing side and of quantifying probabilities in the interesting transition region 
from the dark to the light terrains. As described in Sec. 2.4, wrapping onto 
the trailing hemisphere cannot be the result of lapetus' negligible gravitational 
focusing of retrograde particles (the maximum deflection of a retrograde particle 
by lapetus' gravity is 1°). Such extension is, however, a natural consequence 
of dust particles on eccentric orbits (see Fig. [S]). We therefore now calculate the 
probability distribution in the more general case of eccentric, inclined particles. 
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Figure 9: The orbits capable of striking lapetus are well approximated by a uniform disk of 
parallel trajectories, shown on left. Probabilities are then simply proportional to the projected 
area, given by dA costp. The apex of motion is at the leftmost point on the semicircle. 

Following Greenberg's formalism (1982), we express the probability distri- 
bution function (pdf) as an integral over all the uniformly distributed angles 
{Qi,ujj, Qp,ujp), where the T subscripts refer to lapetus, and the 'p' subscripts 
to the dust particle. Figure 2] shows the geometry of an arbitrary orbit's three 
angular orbital elements i, and lu. 

As Greenberg (1982) notes, however, the problem's geometry does not de- 
pend on the values of Qj and Qp independently — the only geometrically mean- 
ingful quantity is their difference Af2. Furthermore, an important simplification 
can be made by approximating lapetus' orbit as circular (its actual eccentricity 
= 0.03 and for circular dust orbits would lead to extensions of only ^1° onto the 
trailing side). This obviates the need to specify w/, the position of pericenter in 
lapetus' orbit. These two considerations reduce the phase space dimensionality 
from {flj,u!i,ujp,flp) to (wp, Afi), so that 



where the integral spans the region in (ujp, AQ) space in which collisions occur. 
p{9, (j), ujp, AQ) can further be expressed in terms of the conditional pdf p{9, cj)) 
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given the set (wp, Afi) multiplied by the probability density for (cjp, Afl), 

pi9,(i),ujp, An) ^ p{e,,p\Ljp,An)p{LUp, An), (9) 



so 



p{e,(j)) = Jp{e,<j)\ujp,An)p{ujp,An)dujpdAn. (10) 



This simplifies the problem because p{ujp, Af2), the prob ability of strikini 



5 lape- 



19821 ). The 



tus (anywhere) given LOp and AfJ, is already available (|GreenbersJ 
problem is then reduced to finding p(0, (j)\iOp, Ail). 

While analytically correct, Eq. pUj) is a formidable integral to compute 
numerically due to the scale separation in the problem. The orbit is so large 
compared to the satellite that a minute change in LUp or An shifts the location 
of impact drastically. This sensitivity of p(9, (f)\iOp, An) dictates extremely fine 
stepsizes in the integration. When combined with the fact that the calculation 
must be done for 10'* separate timesteps, 8 initial conditions and 6 particle 
sizes, the scale separation indicates a brute force approach will be cumbersome 
at best. 

However, this approach considers each orbital orientation individually. The 
scale separation lets us consider well-defined groups of orientations. Consider an 
orientation where the orbits cross exactly, i.e., the particle would pass through 
the center of lapetus. There is a range in An and ujp around this orbital 
orientation where the orbits no longer exactly cross but are still close enough 
that the particle impacts lapetus (Fig. [TU]) . 

As previously argued, the facts that aj^Rj and that gravitational focusing 
is negligible mean that, near impact, we can approximate these orbits as a 
uniformly distributed disk of parallel trajectories. This approach allows one to 
coarsen stepsizes while retaining the symmetries in the problem and maintaining 
the fidelity of the final distribution. 
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Figure 10: Orbits can only cross along the line that marks the intersection of both orbital 
planes (the line of nodes). lapetus is depicted at one of the nodes, with its size greatly 
exaggerated. For the particle orbits, there is a range in the angle from the node to pericenter 
(ojp) Aojp where collisions with lapetus are possible. Similarly there is a coUisional range in 
AQ, the angle that rotates the line of nodes in the plane (not shown). 



With these considerations, the problem of calculating the distribution func- 
tion is more tractable. Our approach will be to first find the latitude and 
longitude for exactly crossing orbits, and then to find how the disk of parallel 
orbits around the central orbit maps onto the spherical surface of lapetus. 

Given a particular AD,, one can combine it with the inclinations to deter- 
mine the orienta tion of the orbital planes relative to each other (see Fig. 2 in 



Greenberg 



19821 1. The relative inclination i' is given by spherical trigonometry, 

cos i' = cos ii cos ip + sin ij sin ip cos Afl. (11) 

Within the particle's orbital plane, ujp sets the orientation of the orbit. Having 
approximated lapetus' orbit as circular, we do not have to consider the satellite's 
orientation within its orbital plane. From Fig. IIOI it is clear that most particle- 
orbit orientations do not result in a crossing. Furthermore, if the two orbits are 
to cross, they must do so at either of the two nodes where lapetus' orbit pierces 
the mutual line of nodes. 
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The facts that the orbits must cross at a node, and that those respective 
points on the orbits must therefore be equidistant from Saturn sets the pos- 
sible values of ujp (note that this would not be the case if both orbits were 
substantially non-circular). The angle from pericenter to the ascending node 
is, by definition, —uip. From the equation for an ellipse, we therefore have the 
condition. 



ai = 



ap{l 



Rearranging, 



cos U)p 



1 + Cp C0S(— Wp) 



ai 



(12) 



(13) 



Since cosine is an even function, Eq. (jl3l) gives two solutions iwp, reflecting the 
ellipse's symmetry across its long axis. Similarly, another pair of solutions are 
present at the descending node, located at an angle of 180 - iOp from pericenter. 
In genera l, therefore, four orientations yield crossing orbits for a given Af2 (see 



Fig. 1 in 



Wetherill 



1967 ) 



The imposed circularity of lapetus' orbit means that these four orbits will 
strike symmetrically about the equator and about the longitude that corre- 
sponds to the apex of motion (this occurs because the satellite rotates syn- 
chronously). In other words, if we define the longitude in the direction of mo- 
tion as 0°, the four exactly-crossing orbits will strike at (±0, ±0). These can 
later be adjusted to conform with the conventional longitude in the direction of 
motion — the zero-longitude m eridians are currently under revision for Saturn's 



satellites ( Roatsch et al 



120091) . Given a Af2, one therefore only has to compute 
{9, (p) for one of the four orientations and then straightforwardly substitute for 
the other three. Here we choose to consider collisions at the ascending node. 

The relevant vector to consider is the relative velocity vector in a coordinate 
system centered on the ascending node (see Fig. [TT|) . 



The spherical angles that define the direction of the particle's relative ve- 
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Figure 11: lapetus' circular orbit is executed in the lighter horizontal plane, while the par- 
ticle's orbit is carried out in the inclined darker plane, with the two crossing at the particle 
orbit's ascending node. We choose to work in a cylindrical coordinate system centered at the 
ascending node where the z direction is lapetus' orbit normal. For simplicity, we first express 
the particle's velocity in its own orbital plane, where z' is the orbit normal. The relative 
inclination i' and argument of pericenter u)p are also depicted. 



locity vector in this system determine the location of impact on the lapetus 
surface. The relative velocity vector is given by 

Vre/=Vp-V/- (14) 

Having approximated lapetus' orbit as circular, v/ is always azimuthal. v/ can 
therefore be simply written in terms of the uniform circular velocity, 





^.-r^A- (15) 



ai 







Vp is given in the particle's orbital plane in cylindrical coordinates by 



Hamilton 
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(19931), 




Cpsinf 



) 



GMs, 



1 + CpCOsf 



(16) 



ap(l -e2) 







Since we are interested in the particle's velocity at the ascending node in partic- 
ular, we plug in / = —Up (see Fig. [TT|) . At the ascending node, the unit vectors 
r and r' align, but the remaining unit vectors are misaligned by the relative 
inclination i' . We therefore rotate Vp by an angle —i' around the r axis (see 
Fig. [TT]), yielding 



The relative velocity vector in Eq. (ITi|) can then be obtained from Eqs. (|15p 
and pT| . plugging in for cosujp from Eq. p^ . The latitude 9 and longitude 
on lapetus where the particle strikes are then given by. 



The symmetry described earlier can then be used to find the latitudes and lon- 
gitudes for the other three orientations. We have thus determined the location 
where particles on the four possible crossing orbits (for a given Afl) would im- 
pact. The last piece is to include the disk of parallel trajectories around these 
crossing orbits that can still impact lapetus (see Fig. [TO)) . 

As mentioned earlier, the scale separation in the problem allows us to con- 
sider all the orbits that can strike lapetus close to the crossing orbit as parallel 
lines with a uniform probability distribution. The situation is analogous to the 




(17) 



sini'{l + CpCOSLOp) 




(18) 



(19) 
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one presented at the beginning of the section, except with the ip ~ direction 
now interpreted as the incoming trajectory at latitude and longitude 6 and 0, 
respectively. The probability is again proportional to the projected area, so 
normalizing the probability distribution we obtain 

= (20) 

which falls to as V' reaches 90° like it should. Obviously this only applies to the 
hemisphere facing the disk — for wherever ip > 90°, P{tp) — 0. Again, this would 
not be the case with substantial gravitational focusing, but as a result of the high 
relative velocities due to dust orbits being retrograde, gravitational focusing 
is negligible (Sec. 2.4). The probability P{ip) can then be straightforwardly 
converted to a probability per d9 and dcj) for substitution into Eq. (jlOp . 

This provides a prescription for numerically computing p{9, cj)), the probabil- 
ity density function that we originally set out to find, as a function of latitude 
and longitude . Cycling over Ail, at each step, we identify the fo ur crossing 



orbits and their associated probabilities within the interval, using 



Greenberg 



(|l982f ). Then, for each of the four crossing orbits, we "spread" the respective 
probability across the hemisphere defined by the crossing orbit through the 
distribution in Eq. (^0)) . 

Note that, since p{6, (p) depends implicitly on particle eccentricities (cf. Eq. 
[T7|) . it will also be a function of particle size. We can write this explicitly, 
and straightforwardly convert p{9, (p) to a probability per unit area, by defining 
P{r,9,4>) = p{r,9,4>)/cos{9)Ri'^ . We choose to use P{r,9,(j)), the normalized 
probability per unit area, in subsequent calculations. 

3.2. Calculating Depths 

Since radiation pressure produces different orbital histories and is particle- 
size dependent, different particle sizes have different pdfs. Fig. [T^ shows the 
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probability density functions for 5, 10, 50 and 500 /im particles, with each 
contour representing a successive 10-fold decay from the peak value at the apex 
of motion at the left-most point of each figure. 






Figure 12: Moving from the top-left figure clockwise, probability density functions for 5, 10, 
500 and 50 fim particles. Plots represent equatorial views where the vertical line in the center 
represents the boundary between leading and trailing sides. As such, the apex of motion is at 
the leftmost point on each figure. Contours represent successive 10-fold decays from the peak 
value at the apex of motion, down to 10~^ of the apex value. Dot-dashed lines are drawn 
every 10° in longitude. 



The figure shows that smaller particles extend farther onto the trailing side 
near the equator. This is due to their higher eccentricities, as discussed at 
the beginning of Sec. 3 (see Fig. [8]). The distributions for particles > 50/im 
quickly converge to the large-particle limiting distribution, depicted for 500/im 
particles. The eccentricities of these larger particles are too low to cause them 
to significantly wrap around the equator onto the trailing side; however, the 
coverage over the poles is dominated by the precession of lapetus' orbit, which 
is independent of particle size. 

We now use such probability density distributions P{r,9, (f) to estimate the 
depth of dust as a function of position on lapetus' surface. The volume of dust 
particles within dr of size r that lands within an area A on lapetus at latitude 
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and longitude (j) can be expressed as 



Volumeir, 0, </>) = N{r) x P{r, 6*, </>) x A x V{r) 



(21) 



where N{r) is the number of dust particles within dr of radius r generated in the 
outer Saturnian system and V{r) is the volume of a spherical particle of radius 

r. Unfortunately, the current (or past) particle size distribution N(r) is not well 
constrained observationally. We therefore consider a variety of exponents for 
distributions of the form, 



where D is a normalization constant, and ^ is the (negative) power law index of 
the particle size-frequency distribution. Finally, the depth can be estimated (to 
within a packing efficiency factor) as the volume over an area element divided 
by the area of the surface clement. 

The final integration over the range of particle sizes to find the total dust 
depth is complicated by the fact that larger particles, being less affected by 
Poynting-Robertson drag, take longer to reach lapetus from Phoebe. We can 
consider two limiting cases: 

a) Most of the debris in the outer Saturnian system was generated early in the 
Solar System's history (i.e., the mass contribution from the Phoebe ring is neg- 
ligible). In this limiting case, all the particles with lapetus-collision timescales 
(and destruction lifetimes) smaller than the age of the Solar System will have 
had time to impact lapetus and collision timescales across different particle sizes 
are irrelevant. 

b) The mass in the outer Saturnian system has been generated at a constant 
rate over its history. In this case where particles are continuously resupplied, 
smaller particles that decay inward faster will have a larger effect than they 
would have in the first case. 



(22) 
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We begin by considering case a) where we investigate all particles with colli- 
sion timescales tc smaller than the age of the Solar System on an "even footing." 
Since dust particles' semimajor axes decay exponentially through Poynting- 
Robertson drag on a timescale tpr, and since the ratio of Phoebe's to lapetus' 
semimajor axes is « 3.6, particles that strike lapetus do so on roughly a single 
e- folding timescale, i.e., tc ^ Tpn. The P-R timescale is given by Eq. 
This implies that the largest dust size to consider is 1 cm, with corresponding 
TpR ^ 1 Gyr. 

Integrating over all particle sizes, 

Depth{e,(j>) (X f r^-^ X P{r,e,(j))dr. (23) 

For rmirn we use the smallest size of long-lived particles from Phoebe, approx- 
imately 5 /xm (see Sec. 2.1). At the other limit, we use r^ax 1 cm. Should 
lifetimes from cat astrophic collisions between particles or other processes (see 



Burns et al 



20011) be lower than ~ 1 Gyr, rmax must be considered more care- 
fully. 

We now address case b), where particles are continuously resupplied. In 
this circumstance we can consider each particle size to fall onto lapetus at a 
characteristic rate, 

/ „ ,N Volume{r,9,4>) 
i?ate(r, 6*, </>) = L^.^^ ^24) 

where tc is the characteristic collision timescale and is tpr. We can express 
the depth then as 

i?eptM^^,0)oc^^^^^^^;4i^, (25) 

where t is the interval over which dust has been accumulating. From Eq. ([2]), 
we find that tc oc r, so plugging in for the volume as was done in the first case, 
we find that 

Depth{0, 0) oc / r^-'^ x P(r, 6, (j))dr. (26) 
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A comparison between Eqs. (|23p and shows that a constant rate of 
dust production simply acts to steepen the effective power-law index, because 
small particles will arrive at lapetus more quickly than large ones. The effective 
power-law index will be intermediate between the limiting cases of Eqs. (|23p 
and and the depth can therefore be generally expressed as 

Depth{9, 0) oc / r^-iP+^) x P{r, 9, <P)dr, (27) 



where 7 is a number between and 1 that parametrizes the constancy of dust 
production over the age of the Solar System. A 7 of and 1 would therefore, 
respectively, correspond to cases (a) and (b) introduced at the beginning of this 
section (3.2). 

The constant of proportionality in Eq. (|27p is a priori highly uncertain. 
An important, though poorly constrained, quantity is the time by which lape- 
tus had become tidally locked to Saturn. lapetus' dichotomy could not have 
formed prior to this time, as a non-synchronously rotating lapetus would re- 
ceive dust equally on all sides. Furthermore, the thermal models required to 
explain its sharp albed o boundaries and b r ight p o les require lapetus' s l ow 79 



day synchronous period 



Spencer and DenkI pOlOl ) 



Castillo-Rogez et al 



(120071) 



estimate t idal locking occurre d between 200 Myr and 1 Gyr after formation. 



Moreover, 



Bottke et al 



(|2010[ ) argue that most of the dust in the outer Saturn 
system should have been generated in the first few 100 Myr. Case a) reflects a 
situation where most of the dust in the outer Saturn system was generated early 
and lapetus was able to quickly achieve synchronous rotation so that this dust 
mass arrived after locking. If, on the other hand, the timescale for tidal evolu- 
tion is long (^ 1 Gyr), one might expect the production of the relevant dust to 
be fairly constant — case b) — since any initial flurry of dust (should there have 
been one) would have arrived too soon. Despite these uncertainties, one can 
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still normalize the depths over the surface a posteriori through a measurement 
of depth at a particular position {9,(j)). 

In studying lapetus with the Cassini radar instrument, 



Black et al 



t004) 



found little hemis pheric asymmetry in albedo at a wavelength of 1 3 cm, while 



Qstro et al 



Ostroetal 



(120061 ) observed a strong dichotomy at 2 cm wavelength. 
20061 ) interpret these results as implying contamination of ice with dark mate- 
rial to a depth of one to several decimeters. The latter 's measurement on the 
leading side of lapetus was centered on (66° W, +39° N) , but the beam size was 
comparable to the angular size of the satellite (beam/R = 1.36). 

While this measurement is imprecise, it does set the order of magnitude of 
the dark material's depth. Figure [T3] shows depth contours for three different 
choices of effective power law index, /3e// = /? + 7 in (25), assuming a peak depth 
of dust at the apex of motion (extreme left of each figure) of 0.5 m. Figure [HI 
provides depths following the equator and meridian passing through the apex 
of motion at ^ 90° W for the same cases. 

The top graph in Fig. [T3] shows that one should expect extension of dark 
material ^ 20 — 30° onto the trailing side for all expected particle-size distri- 
butions. Only small particles (< 25/im in size), having more eccentric orbits, 
can significantly reach onto the trailing side. As a result, the shallowest effec- 
tive power-law index (/3e// = 3), having fewer small particles, yields a spatial 
distribution that extends onto the trailing side significantly less. 

The bottom graph shows the extension over the poles. Far from Saturn, solar 
torques dominate torques from Saturn's oblateness and cause orbits' angular 
momentum vectors to precess, keeping the inclination roughly constant. Because 
the inclination is set by initial conditions (i.e., Phoebe's orbital inclination), 
and is independent of particle size, the graphs for all three power-law indices 
overlap. The extension over the poles ("latitudes" > ±90° along the meridian 
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jS + r = 5 



Figure 13: Depth contours representing 10-fold decays from the peak value (at the extreme 
left of the figure) for /3e/ / = /3 + 7 = 3, 4 and 5, assuming a peak depth at the apex of motion 
(extreme left of each figure) of 0.5 m. Plots represent equatorial views where the vertical line 
in the center represents the boundary between leading and trailing sides. Dash-dotted lines 
are drawn every 10° in longitude. 

onto the trailing side) is due to both particle-orbit inclinations and lapetus' 
orbital precession. 

One should be careful in distinguishing measured depths of dark material 
(through radar or otherwise) from depths of dust accumulated over lapetus' 
history. If the model of dust deposition and subsequent thermal ice migration 
is correct, the depth of dark material would be the sum of the contributions 
from exogenous dust and from the native lag deposit (see discussion in Sec. 3). 
As long as the depth measured is substantially larger than the expected depth 
of the lag deposit, the distinction is minor. Figure [T51 assumes a peak depth of 
dust of 50 cm-with no impact-gardening, this would imply an actual depth of 
dark material about 20% greater (with a 10 cm lag deposit). Should improved 
measurements of the peak dust depth become available, the contours on these 
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Figure 14: Top graph shows depth vs longitude along the equator for /3e// = /9 + 7 = 3 
(solid), /3e// = 4 (dashed), /3e// = 5 (dash-dotted). Bottom graph shows depth vs latitude 
along the meridian passing through the apex of motion (longitude ~ 90° W). The concept 
of latitude has been extended beyond ±90° along the corresponding meridian on the trailing 
side of lapetus to show the extension of dark material over the poles. 



maps could be straightforwardly rescaled. 

The additional lag-deposit depths arc not included in our modelling; how- 
ever, since exogenous dust acts as the trigger to thermal ice migration, the maps 
above should be good tracers (at low latitudes) of which areas will be dark and 
which will be light. We can therefore attempt to predict the boundary of Cassini 
Regie. 

While the figures show a rapid fall-off in depth on the trailing side, maps 
of lapetus show no gradation in albedo. Areas initially darkened by infalling 
dust and receiving strong insolation become almost completely blackened. As 
a result, the equatorial regions of the leading side appear uniformly dark, while 
dust-free areas of the trailing side and the colder poles, where the ice that 
sublimated at lower latitudes settles, appear about ten times brighter. 

Thus, if one ignores the contours, the plots above argue for a blackened 
leading side extending between twenty to thirty degrees in longitude onto the 
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trailing side for 3 < /3e// < 5. In fact, this holds true for all /3e// > 3 as shown in 
Fig. [121 At low latitudes, this matches maps of lapetus' albedo well(see Figl?]). 
Deposition of exogenous dust therefore neatly explains the boundaries of the 
dark material at low latitudes for the entire range of likely power-law indices 
(accordingly, Spencer & Denk's (2010) model explains the sharp boundaries in 
albedo as well as the bright poles). 



Figure 15: The extension in longitude onto tlie trailing side, chosen as the longitude at which 
the depth falls below 10~^ the peak value, for different values of /3e// = /3 + 7- The discrete 
steps are the result of the resolution of the calculation — 2° in longitude. 



If indeed lapetus was initially darkened by dust from the outer Saturnian 
system, the extension onto the trailing side seems to exclude the shallowest 
power-law indices in the particle size distribution, /? + 7 < 3. While it is 
encouraging that all /3e/ / > 3 are consistent with the observed distribution, the 
small slope in Fig. [12] for f3eff > 3 renders the longitudinal coverage on lapetus 
a comparatively poor indicator of the responsible particle size distribution. 

Fig. [THl shows the distribution for = 3.5, the power-law index for an 



idealized infinite collisional cascade (Burns et al 



2001 ). It assumes a constant 



supply of particles (7 = 1) and artificially accounts for thermal ice migration by 
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brightening the poles down to the observed latitude of ^ ±60° and by completely 
darkening areas with depths greater than 5/im. 

As previously mentioned, 5/im particles are the smallest Phoebe-generated 
particles that would strike lapetus; therefore, this boundary for the anticipated 
depth is roughly where one should expect to transition from uniform darkness 
to the s tochastic dalmatia n patterns observed in the closest Cassini fly by of 



lapetus (jPenk et al 



2010|). This seems consistent with a visual inspection of 



maps (Fig. [7] and see also ([Blackburn et al 



20101) ) and images of lapetus by 



Cassini like the one shown alongside in Fig. [1^1 




Figure 16: Model of the lapetus surface assuming a peak depth of 50 cm, /3 = 3.5 and 7 = 1. 
All areas with depths > 5fim have been uniformly darkened and the poles beyond ±60° 
latitude brightened to artificially account for thermal ice migration. On the right is an image 
taken by Cassini at roughly the same orientation for visual comparison (obtained from the 
Planetary Photojournal— PIA08273) . 



An important difference between the modeled surface and the observed dis- 
tribution is that the theoretically derived dark terrain is concave, in the sense 
that if one follows a meridian on the boundary, one sees that the dark material 
extends further in longitude at higher latitudes. Cassini Regio is convex, as 
can be clearly seen in Fig. [71 which uses a simple cylindrical projection where 
meridians appear as vertical lines. Since temperatures should drop smoothly 
as one moves from equator to pole, perha ps the discrepancy result s from ther- 
mal ice migration. Indeed, such models 



Spencer and Denk 



20101 ) are able to 
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reproduce this concavity. 



Denketal 



Apart from the long-known albedo dichotomy on lapetus, 
recently detected a new color dichotomy on lapetus in which the leading side of 
lapetus is substantially redder than the trailing side. The color dichotomy seems 
to extend farthe r pole-ward, and t ransitions more gradually onto the trailing 



side. Perhaps, as 



Denketal 



(120101) point out, the color dichotomy more faith- 
fully traces where the dust landed while the albedo dichotomy reflects thermal 
ice migration's modification of the initial pattern of dust deposition. Further 
work is needed to ascertain quantitative agreement between observations of the 
color dichotomy and theoretical models like those presented here. 

3.3. How Much Sublimation Is Required to Paint the Poles Bright? 

Substantial amounts of dust should have struck lapetus at high latitudes; 
however, the poles appear bright (Fig. [7]). If the preceding section is correct, 
this means that bright, sublimed ice from lower latitudes must be settling on 
the polar regions faster than dark dust is landing on them. Our coUisional flux 
then provides an opportunity to constrain sublimation rates. 

We found in the previous section that the distribution at the boundary 
between the leading and trailing sides depends on the underlying particle size 
distribution of dust; however, farther from the boundaries, as argued at the 
beginning of Sec. 3.1, the depth (at low latitudes) should scale approximately 
as cos^ independent of the particle size distribution, where V' is the angular 
distance from the apex (cf. Eq. ((201) '). 

Therefore, following the meridian that passes through the apex of motion 
(longitude w 90°), at a latitude of 60° the depth should be roughly half that at 
the apex (cos60° = 1/2). This point in the polar region would receive more dust 
than any other point at the same latitude as it has the minimum angular dis- 
tance from the apex. The fact that this point in the polar region with maximum 
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dust flux appears bright provides the strongest constraint on the minimum sub- 
hmation rate required to keep the poles bright. Assuming a peak depth of dust 
at the apex of 50 cm as done in the previous section based on radar measure- 



ments ( Ostro et al 



I2006I ). this implies an average rate for the minimum polar 
dust deposition of ~ 25cm /5Gyr or ^ 50fim/Myr. 

It is possible, however, that the average rate of deposition would not match 
the rate of sublimation. Maybe in the past (w hen deposition rate s of dust in 



Bottke et al 



2010l ). the poles 



the outer Saturnian system were likely higher, 
of lapetus were dark. Perhaps only recently did deposition of ice exceed that 
of dust and hide evidence of past dark poles. In that case the rate of ice 
sublimation required to keep the poles bright would be lower than the average 
dust deposition rate. Dark-ringed craters in the polar terrains could support 
such a conjecture, but current observations are unable to distinguish between 
these two general possibilities. 

Therefore, given only the observation that lapetus' poles are bright today, we 
now try to roughly constrain the current sublimation rate. We can estimate the 
deposition rate of the material coming from the Phoebe ring using its measured 
optical depth. 

For lack of better information, we assume that the entire volume of the ring 
has the same particle size distribution. In this case, the rate of dust deposition 
at a latitude of 60° along the meridian passing through the apex (longitude 
- 90° W) can be directly obtained from Eqs. (p4)) . ((2T|) . ([22]) and ([2|). 



Dr-i^ VdP{r, 60°, 90°) ^ttt^ 
0.1 r ^ 



Rate{m° ,90°) ^ / ' ' ^ dr, (28) 



where Vd is the volume of the disk and P{r, 60°, 90°) is the probability per unit 
area for dust striking at the longitude corresponding to the apex (^ 90°!^) 
and a latitude of 60°. In pursuing an order of magnitude estimate, we take 
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all particle sizes in the Phoebe ring (i.e., > 5/ini) to have probability ^ 1 of 
striking lapetus (cf. Fig. [S]). Furthermore, we consider the depth distribution 
to be simply proportional to costp (Fig. [U, with depths of on the trailing side. 
For the integrated probability to yield unity, P(r, 60°,90°) w 3 x 10^^ km^^. 

We esti mate Vd as the volum e of a disk 300 Rs in radius and 40 Rs thick 
(as done by 



Verbiscer et al 



20091) . or - 2 X lO^^ km^. This yields 



Rate{60°,90°) - 10 



km, fim f^" 



Myr 



(29) 



The final integral here is provided by the definition of the normal optical 
depth. 



n{r,l)QextCT{r)drdl, 



(30) 



Verbiscer et al 



(|2009| ). we take values 



where Qext = Qabs + Qscat- Following 
of Qabs — 0.8 and Qscat = 0.2. They estimate t 2 x 10~^. Making again 
the simplifying assumption that the number density does not depend on the 
distance along the line of sight / (trivializing the integration over I), 

2 X IQ-s 



Dr^-Pdr. 



(31) 



40i?S TT 

If we are interested in the rate of deposition due to the material currently 
seen in the Phoebe ring (i.e., if we take the same limits of integration), the 
integral can then be plugged into Eq. (|29p , yielding a rate of dust deposition at 
the leading edge of the polar region of ^ 50fi'm/Myr, which is of the same order 
as the averag e deposition rate calculated earlier. This estimate also agrees with 



that given by 



Verbiscer et al 



(120091) of ~ AOfim / Myr. 
In either event, lapetus must today be actively depositing on the order of 
tens of fim per Myr of sublimed ice onto the poles in order to keep them bright. 



This ice could originate from either the leading or trailing side. In Cassini 
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Regio, with daytime temperatures of 130 K, i ce will sublime over ext remely 



20101) . Very 



short timescales — about 1000 /im in 8000 years (jSpencer and DenU I; 
quickly, a dark layer thicker than the thermal skin depth will form, making it 
difficult for further ice to sublimate. In this case, the rate of impact gardening, 
which brings fresh ice to the surface, will determine the rate of sublimation from 
Cassini Regi o. Unfortunately, impact gardening depths on lapetus are not well 



constrained ([Spencer and Denk 



2010). 



The temperatur e on the brighter trailing side indicates a sublimation rate of 



about IQQfjLm/Myr (jSpencer and Denk 



20101 ) . Given that the area of the trailing 



side at latitudes < 60° is comparable to that of both polar regions combined, this 
sublimation rate seems capable of overwhelming polar dust deposition. Further 
data on the Phoebe ring will help further constrain the required sublimation 
rates, and improved modeling of the surface processes on lapetus will ultimately 
dictate the consistency of these two pictures. 

4. Dust from other Irregular satellites 

A long-standing objection to Soter's model of dust infall from Phoebe has 
been that Cassini Regio's spec t rum differs from tha t of Phoebe (e.g. 



Tholen and Zellner 



1983 : 



Cruikshank et al 



1983 



Buratti et al 



20021 ). One possible resolution is 



that the surface co mpositions are in fact similar, but other effects such as 



Rayleigh scattering (jClark et al 



20081 ). cause the spectra to differ. Another 



possibility is that lapet us has also been coated by dust from the irregular satel 



lites other than Phoebe (IBuratti et al 
and 



Buratti et al 



2005 



Tosi et al. 



2010) 



Grav et al 



mm 



(j2005l ) find that many of these new irregular satellites have a 
reddish color similar enough to that of Hyperion and Cassini Regio to suggest 
a link between them. 

Unlike the regular satellites (ignoring lapetus) that move on low-eccentricity 
orbits close to Saturn's equatorial plane, the irregulars have widely varying 
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inclinations and eccentricities. This implies a violent collisional history between 
irregular satellites as differ ing precession rates w ould have led to crossing orbits 



and consequent collisions (iNesvornv et al 



collisional grinding numerically, 



Bottke et al 



2003). B y modeling this process of 



(|2010D estimate that - lO^o kg of 



dust should have been generated in the outer Saturn system, particularly early 
in the Solar System's lifetime. 

We therefore explore the likelihood that debris from these other irregular 
satellites would collide with lapetus. Fig. [TTl shows the probability that 10/im 
grains will strike lapetus if they start with the orbits of the various irregular 
satellites known today, plotted against both today's value of the parent-satellite 
orbit's inclination and eccentricity. The orbits of the current irregular satellites 
are not chosen to be necessarily representative of their orbits over the course 
of the Solar System's history — the irregular satellites seen today are likely the 
fragments of past satellites and are subject to increasingly strong gravitational 



Nesvornv et al. 


2003 


Turrini et al. 


2008 


Bottke et al. 


2010) 



2010 ). Rather, we chose 



the current orbits as a way of sampling the orbital phase space of irregular 
satellites. 

The efficiency with which material is supplied to lapetus differs markedly 
between the prograde (plus signs) and retrograde (diamond) satellites (Fig. 
I17p . This can be understood through a simple particle-in-a-box estimate of 
the collision timescale, where the irregular-satellite and dust-particle orbits are 
taken to precess around the same axis (normal to Saturn's orbital plane) : 



Tcoi ^ sin^ip + sin^ij 



ai 
Ri 



Ur 



(32) 



where the p subscript refers to the dust particle, ip and ij are inclinations 
measured relative to Saturn's orbital plane, Rj is lapetus' radius and a/ its 
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Figure 17: Collision probability with lapetus for lO/im grains that start with the orbits of 
today's irregular satellites. Plus signs represent prograde irregulars, while open diamonds 
are retrograde. Inclinations are measured relative to Saturn's orbital plane. The asterisk 
represents Phoebe, and the boxed diamond Ymir (of importance below). 



semimajor axis; U is the relative speed between the two objects, and Ur is 
the radial com ponent of the relative velocity; finally, To^b is the dust particle's 



orbital period (jOpik 



1951 



Hamilton and Burns 



1994I ). lapetus' orbit does not 



quite precess around the normal vector to Saturn's orbital plane, but rather 
around an a xis (the normal to its local Laplace plane, see Fig. [5]) ~ 11° away 



(Ward 



198lh : however, for our rough estimate this can be ignored. 
Qualitatively, one should expect Tcoi to decrease as ip approaches 0° (or 
180°). As ip approaches coplanarity, the phase space that the particle must 
explore before "finding" lapetus is reduced. One should also expect prograde 
particles to have a substantially decreased chance of striking lapetus compared 
to retrograde particles, as retrograde particles have a much larger relative veloc- 
ity U than prograde particles. This occurs because for larger relative velocities, 
when one particle is passing through the node where collisions are possible, the 
other particle can initially be at a wider range of positions in its orbit and still 
reach the node "in time." 
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For prograde particles with low inclinations, the azimuthal component of U 
is mostly subtracted o ut so Ur/U is rouRhly 1, a nd nearly independent of the 



particle eccentricities (jHamilton and Burns 



1994I ). For retrograde particles, on 



the other hand, the azimuthal component of U dominates so Ur/U will be small 
and will increase with orbital eccentricity, which determines the departure of the 
dust orbit from being purely azimuthal. One should therefore expect that for 
retrograde particles, those with inclinations closest to coplanarity and with low 
eccentricities will have the shortest collision timescale and the highest collision 
probability (see Fig. [T7)) . 

We can also investigate this more quantitatively. The ratio Ur/U can be 
obtained from Eqs. and (15), yielding 



Ur^ 

IT 



^ ^ a2(l + ^2 _ 2a cosi') 



(33) 



e2 _ (q,2 _ 1)2 

where i' is the mutual inclination between the particle's and lapetus' orbits, e 
is the particle's eccentricity, and 

a' = (34) 

ai 

We note that while particle orbits will perform small oscillations in their in- 
clinations around their parent body's inclination, particle eccentricities will be 
substantially larger than those of the parent bodies due to radiation pressure. 
In this section we therefore take ip ~ isateinte and select characteristic eccen- 
tricities from our numerical integrations. It is nevertheless generally true that 
source satellites with more eccentric orbits will yield particle orbits with higher 
eccentricities. 

Since particles can only collide with lapetus when their orbit's pericenter is 
smaller than aj and their apocenter is greater than a/, ranges between 1 — e 
and 1 + e. Furthermore, from Eq. i' varies between ip — ii and ip + ij. Taking 
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i' K, ip and a ~ 1 as characteristic values, and expanding cosi' to leading order, 

/ \ o 



Ur 
IT 



(35) 



Since sin^ij <iC siri^ip for Saturn's prograde irregular satellites (35° < ip < 
50°), we can approximate the first term in parentheses in Eq. (1321) as simply 



sinip K, ip. Therefore, 



col 



{prograde) . 



(36) 



Over the range of prograde inclinations and the range of characteristic particle 
eccentricities found in our numerical integrations (0.3 ^ e < 0.6), Tcoi varies 
by a factor between ^ 0.3 — 0.45. This approximation agrees with the values 
from Eq. (prograde) to within 25% over the range of characteristic values for 
e, ip and Op for dust from the prograde satellites, and matches the small spread 
(< 50%) in collision probability for the prograde satellites (plus signs) in Fig. 



For the retrograde satellites, again taking i' « ip and awl, and noting that 



TT — is small so that cosip 



-l + (^- V)V2, 



T, 



col 



4(1-^^) 



(37) 



Thus, since (tt — ip)^/4 <C 1 for the retrograde satellites and 4/e^ ^ 1 for the 
range of characteristic eccentricities (> 0.3), 



Tcoi ~ (tt - ip) - {retrograde). 



(38) 

This agrees with the values from Eq. (pS)) to within 20% over the range of 
characteristic values of e, ip and Op for dust from the retrograde satellites. It 
also shows why the retrogrades have much higher collision probabilities. For the 
range of dust inclinations and eccentricities, the prograde to retrograde ratio of 
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Tcol using Eqs. ([5^ and ([55]) is always greater than unity and < 25. 

Eq. ([55)) means that low-eccentricity moons (those that yield lower eccen- 
tricity particles) with inclinations close to 180° will have the shortest collision 
timescales, and therefore the largest collision probabilities with lapetus. Fig. [T5] 
shows the same probabilities as Fig. 1171 but in two dimensions so as to separate 
the dependence on inclination and eccentricity. Darker squares represent larger 
collision probabilities. Thus one can see following rows of constant eccentric- 
ity that the probability increases as the inclination approaches 180°, whereas 
following colum ns of constant incli nation, the probability decreases with increas- 
ing eccentricity. 



Tosietal 



(|2010l ) find similar trends using a different method 
of evaluating collision probabilities. The two low-eccentricity, high-inclination 
moons on the bottom right of the plot with the highest probabilities are Sut- 
tungr and S/2007.S2. 

4.1. Dust Generation Efficiencies 

While the previous section addressed the likelihood of particles from differ- 
ent irregulars striking lapetus once they are ejected, one must still determine 
the relative dust yield from the various satellites to infer the dominant sources 
of dust for lapetus. Dust will be generated both in collisional break-up between 



the outer irregular satellites (jBottke et al 



2010) and in micro meteoroid bom 



19991) . In both cases. 



bardment from outside the Saturnian system (jBurns et al. 
the effectiveness of a satellite as a dust source is determined by the competition 
between a larger satellite radius raising the collision cross section and a larger 
satellite mass increasing the escape velocity (thus inhibiting dust from leaving 
th e satellite). 



Burns et al 



(|l999r ) investigate this relationship in the Jovian system. For 
small moons, where gravity is not important, the rate at which mass is supplied 
to the ring by a satellite of radius Ri is simply proportional to Ri^\ however, 



48 



1 



;f 0.75 

c; 
u 
u 

LU 

2 0.5 

la 
L_ 

O 



0.25 



130° 140° 150° 160° 170° 180= 
Satellite Orbital inclination 



Figure 18: Numerically computed probabilities for lO-^im dust particles striking lapetus as 
a function of parent (retrograde) satellite orbital inclination and eccentricity. Probabilities 
range from darkest (Suttungr = 0.9 and S/2007_S2 = 0.89) to lightest (Narvi = 0.22). Collision 
probability increases as the inclination approaches coplanarity (180°) and as the eccentricity 
decreases. 



beyond an optimum satellite size Ropt that depends on regolith pr operties, the 



dust 



Droduction rate becomes almost flat, decreasing as Ri 



Burns et al 



(|l999f ) estimate that this optimum size should be about 5 — 10 km in the Jo- 



vian system. Assuming similar results for the Saturn system, this implies that 
Phoebe [R ^ 100 km) should produce no more (in fact slightly less) dust than 
any ~ 5 — 10 km irregular satellite. Fig. [TH] shows the impact probabilities 
calculated in the previous section for each irregular satellit e weighted by the R j 



Burns et al. 



mm 



dependent terms in the dust-generation efficiency factor of 
assuming an optimum satellite size Ropt of 10 km. 

Fig. [19] shows that Ymir {Ri = 9 km) should be roughly as important a 
contributor of dust to lapetus as Phoebe {Ri ~ 107 km), though the summed 
contribution from the remaining moons is greater than that of either Ymir or 
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Figure 19: Dust supplied to lapetus from each of the irregular satellites relative to the con- 
tribution from Phoebe, plotted vs. satellite radius. The relative contribution is calculated as 
the product of the collision probability for the particular satellite (Sec . 4) a nd the radius- 
dependent terms in the dust-generation efficiency factor of iBurns et al.l lll999l ). assuming an 
optimum satellite size of 10 km. Prograde satellites are represented by plus signs, and retro- 
grade moons by open diamonds. Ymir is plotted as a boxed diamond {Ri =9 km) and has 
the largest contribution. 



Phoebe. This might help lessen the contradictio n that the spectrum of Cassini 



Regio does not seem to match that of Phoebe (jBuratti et al 



2005 



Tosi et al 



2010h . 

But if Phoebe is not the dominant source of dust in the outer Saturnian 
system, why then is the only prominent dust ring generated by the irregu- 
lar satellites associated with Phoebe? Satellites should generate dust rings of 
height 2asin i, so one might expect to see a nested series of rings of differing 
heights. This is analogous to the dust bands observed in the zodiacal cloud, 
where the dimension s of the bands give aw ay the orbital elements of the object 



that produced them (jPermott et al 



1984 ) 



Perhaps the fact that Phoebe's orbit has the smallest semimajor axis and 
lowest inclination among the irregular satellites squeezes its modest share of 
dust into a more compact volume, yielding a higher optical depth than other 
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satellite rings. The line-of-sight optical depth of a ring generated by satellite j, 
Tj UjajLj where rij is the number density, aj is the average particle cross- 
section and Lj is the distance along the line of sight. Focusing only on the 
parameters involving the satellites (as opposed to dust properties), and taking 
Lj aj, Tj ~ MjUj/Vj, where Vj is the volume of the ring associated with 
satellite j and Mj is the mass of dust within it. Vj should be proportional to 
a^sinij, yielding 



Mi 



ajSm lj 



(39) 



The mass Mj carries the same weighting factor used in Fig. [121 The optical 
depth therefore carries an additional factor of {ojSinij ) . Fig. [201 plots the 
weighting factor in Eq. p9p relative to that for Phoebe vs. the expected ring 
height that would be produced by the particular moon {2asini) in Saturn radii. 
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Figure 20: Estimated line-of-sight optical depth of rings created by the irregular satellites 
relative to the optical depth of the Phoebe ring, plotted vs. the height of the ring that the 
satellite would produce in Saturn radii. Phoebe (asterisk) generates the ring of highest optical 
depth with a thickness of ~ 40 Rs, followed by Ymir (boxed diamond), which should produce 
a ~ llOiJs-tall ring. Plus signs denote regular satellites, open diamonds irregular moons. 



The simple estimates illustrated in Figs. [19] and [20| suggest that irregular 
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satellites other than Phoebe could contribute substantial amounts of dust to 
lapetus, while the Phoebe ring (due to its compactness) would be the most 
prominent dust ring generated by the irregulars. Perhaps in observations of the 
Phoebe ring, the flux from the various much taller, low-optical depth rings has 
been interpreted as part of the background. One might still, however, expect 
to be able to identify a dust ring associated with Ymir; it would be 3 times 
taller, and have ~ 1/4 the line-of-sight optical depth of Phoebe's structure. 

The fact that no other dust rings have yet been detected could be used to 
argue that other factors enhance Phoebe's dust production in the outer Sat- 
urn system. One possibility is that Phoebe's position as the innermost irregular 
satellite increases its c ollision frequency with other irregulars. Numerical studies 



( Nesvornv et al 



20031) suggest that Phoebe alone among the Saturnian irregu- 
lars likely suffered collisions with several now absent irregulars; while below the 
detection threshold of today's telescopes, the ejecta from these events would be 
excellent suppliers of debris. As such, there might be an increased amount of 
unseen coUisional debris (< 1 km) sharing Phoebe's orbital elements, which, for 
a steep enough size distribution, could contribute significantly to the Phoebe 
ring. A more certain assessment will have to await further observations of the 
Phoebe ring and searches for separate dust bands. Our studies assume that 
Phoebe is the dominant dust source in the outer Saturn system. Should evi- 
dence to the contrary arise, further work would be required to assess the relative 
contributions to lapetus, Hyperion and Titan. 

5. Implications Beyond lapetus 

5. 1 . lapetus as a Tracer of the Initial Dust Mass at Saturn 

(|2010r ) argue that most of their estimated 10^° kg of collisionally 



Bottke et al 



generated dust in the outer Saturn system was created within a few hundred Myr 
of the capture of the irregular satellites. In this case, the whole range of particle 
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sizes we considered should have had time to reach lapetus, corresponding to the 
7 ~ case discussed in Sec. 3.2. In this circumstance, since particles roughly 
larger than lO/im are almost certain to strike lapetus, we should expect all the 
dust mass in sizes > 10/zm to be p a rt of Cassini Regio. 



Unfortunately, as 



Bottkeetal 



pOlOl ) point out, 10^*^ kg would generate a 
dark layer on lapetus that is kilometers thick. Taking the depth on the leading 
side to fall off as cos ip, we can estimate the volume of dust on lapetus as 

dn 



Via 



TtR^ da 



850 



^apex 

50cm 



(40) 



where dapex is the peak depth of dust at the apex (not including any lag deposit 
from sublimation). Assuming 100% transfer efficiency (all particles > lO/im), 
this would imply an initial dust mass ^ 10^'"' kg. 

The transfer efficiency could be reduced with a sufficiently steep size dis- 
tribution (one that would cause essentially all of the dust mass to be in sizes 
< 5/xm). These small particles would then be quickly eliminated from the sys- 
tem through radiation pressure and avoid lapetus; unfortunately, the requisite 
power-law index is implausibly steep (> 5.5). 

The two results might also be reconciled if lapetus achieved synchronous 
rotation much later (> IGyr after formation). In this case, the dark material 
would have been localized to the leading side only after the bulk of the influx had 
occurred. Such a scenario could help explain why even the bright hemisphere 
of lapetus is darker than the surfaces of the other large icy satellites; however, 
it also poses the problem of why the masses of dust that would have blanketed 
lapetus at all longitudes prior to sync hronous rotation did not cause blackening 



everywhere when thermal migration (jSpencer and Denk 



20101 ) kicked in. 



It seem s for now that, un less the depth estimate derived from radar mea- 



surements (Ostroetal 



20061 ) is grossly in error, the amount of dark material 
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on lapctus implicates an initial dust mas s in the outer Saturn ian system about 
5 orders of magnitude lower than that of [ 



Bottke et al 



POIOO . While the pres- 



ence of ammonia could make a thicker layer of dark material appear shallow 



( Ostro et al 



20061 ). the discovery of small bright- floored craters close to the 



boundary of Cassini Regio support the idea of a thin dark deposit (jPenk et al 



2010 ) 



5. 2. Hyperion 

During Cassini's close fly-by, 



Cruikshank et al 



mm and 



Thomas et al 



(|2007r ) found Hyperion to be segregated into a low-albedo unit mostly filling 
the bottoms of cup-like craters and a more widespread high-albedo unit. Fur- 
thermore, the spectra of the dark material show sim ilarities to the material 



making up Cassini Regio, suggesting a common source (|Cruikshank et al 



Buratti et al. 



2005) 



Burns et al 



2007 



19961) . in their dynamical study of the fate of Phoebe dust, had 
already argued that Hyperion should receive a significant fraction of Phoebe dust 
grains. In their calculations, however, they considered collisions with lapetus, 
Hyperion and Titan sequentially, when in fact — because of radiation-pressure- 
induced orbital eccentricities — dust grains can reach all three moons nearly si- 
multaneously. As a result, Hyperion will receive a much-reduced fr action of the 
grains (~ 0.004 for both 5- and 10-/im grains vs. 0.18 according to 



Burns et al 



(|l996( )). Furthermore, no grains > lO/xm reach Hyperion. The conclusion, how- 
ever, is the same — Hyperion's surface layers should contain some dark material 
from Phoebe and the other irregular satellites. Furthermore, since Hyperion is 
chaotically r otating rather tha n tidally locked, an isotropic distribution of dust 



is expected (jBurns et al 



19961) . 

We estimate the volume of material striking Hyperion relative to the volume 
hitting lapetus as 
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V, 



Hyp Jbfim 



lap 



J IQpm 



(41) 



where PHif) is the probabihty for a particle of size r striking Hyperion (derived 
from our numerical simulations discussed in Sec. 2), Viap is given in Eq. (|40p. 
and we have approximated the probability of particles striking lapetus as a step 
function at r = 10/im. 

Calculated dust depths on Hyperion for 2 < /? < 4 in the limiting cases of 
7 = and 1 (cf. Sec. 3.2) are given in Fig. W\\ assuming isotropic coating and 
a spherical target of radius 135 km. 

Average Dust Depth in mm* 
10 




4.0' 



Figure 21: Average dust depth on Hyperion (in mm) vs power-law index /3, for the limiting 
cases of 7 = (bottom curve) and 1 (top curve) . 



These extremely shallow depths of < 1 cm render it plausible that Hyperion 
might not be uniformly covered in dust, though the mechanism for segregating 
dark material to the bottom s of Hyperion's ubiquito us sharp-edged craters as 



is observed remains unclear ( Cruikshank et al 



20071) . Constraints on the par- 



ticle size distribution in the Phoebe ring from future observations may narrow 
estimates of the material delivered to Hyperion. 
5. 3. Titan 

As discussed in Sec. 2.4, Titan's much larger cross-section causes it to ef- 
ficiently sweep up almost everything that crosses its orbit. The slow inward 
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migration of dust particles > lO/im gives lapetus enough time to collect most 
of them before they become Titan-crossing; our numerical simulations of Sec. 
2, however, indicate that ^ 70% of 5/im particles strike Titan (and that smaller 
particles are so affected by radiation pressure that they either strike Saturn 
in the first half-Saturn year or escape the system). As in the case with lape- 
tus, particles will strike Titan on its leading side; however, its atmosphere will 
distribute material around the entire moon and will fragment particles upon 
entry. 

Using Cassini's measurements, 



Porco et al 



(|2005l ) report a detached haze 



layer at 500-km altitude on Titan. Since the sedimentation time in this layer is 
short, some pro cess must continually repl enish particles. Two hypotheses have 



been proposed (jTomasko and West 



2003) : the la yer either represe nts a con 



densation region at a local temperature minimum ( Liang et al 



2007) or occurs 



20091) . 



where aerosols produced at higher altitudes settle (|Lavvas et al. 

While this haze layer is likely due to the above-mentioned atmospheric ef- 
fects, we can explore whether exogenous dust-infall might also be a significant 
contributor. We therefore derive the volume of particles striking Titan as we 
did for Hyperion (Eq. I41[) assuming a present rate of dust production that is 
constant in time (7 — 1). Fig. shows the results, expressed as a mass flux 
computed ~ 500 km above the Titan surface (i.e., R = 3100 km). 

The cal culated mass flux falls well short of the estimated 2.7-4.6 xlO^^^ 
gcm~^s~^ (jLavvas et al. 



120091 ) required to replenish particles. Fortunately, the 
mechanisms listed above seem sufficient for explaining the haze layer (R. A. 
West, personal communication, 2010). 



6. Conclusion 



Our results show that out of the dust particles collisionally generated at 
Phoebe that are long-lived (grains > 5/im), most larger than ^ lO/im will strike 
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Figure 22: Dust flux into the Titan atmosphere at 500 km altitude (in g cm ^ s ^) vs. 
power-law index /3 for the limiting case of 7 = 1. 



lapetus due to modifications of their orbits by perturbations from the Sun. The 
latter include Poynting-Robertson drag, solar radiation pressure and the Sun's 
tidal gravity in the Saturn system. 

Our computed dust coverage on the lapetus surface matches up well with 
the newly discove red color dichotomy on lapetus that extends up to the poles 



(Denketal 



20101 ). The calculated distribution also traces the shape of Cassini 
Regio well in the longitudinal direction, but realistic thermal modeling is re- 
quired to explain both the bright poles and the sharp boundaries between bright 
and dark material. 

Our orbital calculations for 10/im particles show that dust launched from 
other retrograde outer irregular satellites can have comparable likelihoods of 
striking lapetus to those of dust launched from Phoebe. We argue this can 
contribute to the differing spectra between Phoebe and lapetus; however, the 
question of how much dust was generated by Phoebe relative to the other irreg- 
ular satellites is still unclear. 

By tracking the dust that strikes Hyperion, we find that just a veneer should 
have been laid down on its surface (< 1 cm on average). This picture may 
be consistent with the observation that the surface is not uniformly coated. 
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Cruikshank et alj ( 2007 ) and I Thomas et alj ( 20071) find the dark material to be 



predominantly at the bottoms of Hyperion's ubiquitous cup-like craters. As 
opposed to lapetus, which is tidally locked, Hyperion rotates chaotically, which 
can explain the presence of dark material throughout the surface. 

We determine that effectively all long-lived dust particles that avoid lapetus 
(i.e., a fraction of those between 5 and 10/im) are swept up by Titan. While 
these constitute a considerable mass flux into the Titan atmosphere, they are 
insufficient to account for the satellite's detached haze layer. 
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